1 



Running Coupling in SU(3) Yang- Mills Theory 
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We report about our ongoing computation of running coupling constants in asymptotically free theories using 
the recursive finite size scaling technique. The latest results for the SU(3) Yang-Mills theory are presented. 



1. INTRODUCTION 

In nonperturbative evaluations of asymptoti- 
cally free theories all input parameters can be 
fixed by matching low energy physical quanti- 
ties, like the mass spectrum in QCD, with ex- 
periment. A physical renormalized coupling con- 
stant at high energy is then a computable num- 
ber. If such a computation is achieved with suf- 
ficient precision for some suitable coupling a{q), 
then the connection with the whole perturbative 
sector of the theory can be made by expanding 
in powers of a{q). One has then established a 
link between the high and low energy sectors. Of 
course, with only a few perturbative orders avail- 
able, it is essential to use a non-pathological a{q) 
and to control its values over an adequate range of 
high energies. It is important to understand that 
such a computation of absolute numbers for a{q) 
— even at large q where it evolves perturbatively 
with q — is a nonperturbative problem. As we 
fix parameters at low energy, we have to specify 
g as a dimensionless multiple of some low energy 
quantity like a mass m or the string tension in 
the quenched theory. The difficulty of such a cal- 
culation is mainly due to the fact that the ratio 
of scales q/m should be large, at least 0(10). In 
addition, at least with conventional approaches 
(see below), the lattice spacing o and system size 
L have to be remote of either m or q. 

A straightforward idea to compute the coupling 
in pure gauge theory is to proceed via the static 
quark-antiquark force F(r). While it saturates 
to the string tension (the only free parameter) at 
large distance, F[r) can also be used as a physi- 
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cal running coupling constant at small separation 
r. This requires control of the force from short to 
long distance in one simulation with a <^ r <^ L 
holding for the whole range of physical r involved. 
With L/a always limited to feasible lattice sizes 
like 32 or 48, compromises on the above condi- 
tions have to be accepted, and it is hardly possible 
to vary all scale ratios significantly to check for 
the stability of the results. State of the art calcu- 
lations along these lines are reported in [1,2]. It 
has to be noted that the highest physical energies 
that can be reached here are below about 
2 GeV, if one only stays a factor 2 ... 3 away 
from the cutoff energy. Cutoff effects are cor- 
rected semi-empirically using the lattice Coulomb 
propagator. While these are difficult and careful 
simulations, we find it somewhat hard to assess 
the systematic errors in a completely convincing 
fashion. 

An alternative attempt to derive the coupling 
in QCD has been pioneered by the Fermilab group 
[3]. Here, in a quenched simulation, the spin av- 
eraged IP- IS charmonium splitting is determined 
on a physically large lattice. Although this is a 
nice experimentally known scale with little sen- 
sitivity to the quark masses, also other masses 
could in principle be used here to set the scale. 
The point relevant in the present context is, that 
they extract from such a simulation the bare lat- 
tice coupling go together with the corresponding 
lattice spacing o in GeV. A perturbative method 
is then used to relate go to a physical coupling 
at a scale of the order of the cutoff. The scale 
problem is clearly alleviated in comparison to the 
quark force method, as effectively the cutoff 
is identified with the high energy physical scale. 
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The problem is now to convert reliably from the 
non-universal unphysical lattice-coupling go to a 
physical continuum one like as(q = tt/o) with 
one non-trivial perturbative expansion coefficient 
(presently) known. In [4] it is shown that a large 
number of numerically known quantities, Wilson 
loops for instance, are well approximated by the 
so-called improved perturbation theory. This is 
used in [3], with the average plaquette as the key 
input parameter to implement improved pertur- 
bation theory. Again, one may feel that there 
could be a problem with estimating the error in 
the conversion. 

In the following section I recall the recursive 
finite size method and present the latest results 
obtained. 

2. RUNNING COUPLING BY RECUR- 
SIVE FINITE SIZE SCALING 

The ideas and techniques behind this approach 
were presented in detail in Liischer's talk in Am- 
sterdam [5], and we content ourselves here with 
a brief reminder. The freedom which coupling 
to compute is exploited by choosing a{q) = 
g^{L)/4:'!r,q = 1/L, a coupling that runs with the 
boxsize L. It's definition is based on an abelian 
background field that is induced through nontriv- 
ial boundary conditions [6]. A physical change in 
the free energy with a variation of the field leads 
to the coupling g^[L). The definition for general 
L is independent of perturbation theory, and at 
small L we expect 3 to be a smooth perturbative 
coupling. As in the Fermilab approach one of 
the scales unavoidable in numerical treatments is 
used for physics here to eliminate one large scale 
ratio. We think, however, that L has the advan- 
tage over o of finite size effects being universal. 

A further characteristic of the method is to 
break the problem into several steps. In each step 
of a series of simulations we numerically answer 
the question: Given g^{L) = u, with u some num- 
ber like 1.234 implicitly defining the scale L, what 
is g^{2L)? This result equals the value of the step 
scaling function. 



(t(2,w) = lim i;(2,w, o/L), 



(1) 



izations with decreasing o that are to be extrapo- 
lated to the universal continuum function cr. The 
values, for which cr is determined are tuned such 
that we can iterate cr(2,cr(2, ...)) and thus recur- 
sively get cr for scale factors larger than two. In 
each individual step we only have to deal with 
scales o and L, and we use all the available L/a 
to determine and extrapolate the cutoff" depen- 
dence. All our results for S in the pure SU(3) 
gauge theory [7] are shown in Fig. 1. 
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where S refers to a sequence of finite lattice real- 
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Figure 1. Extrapolation of the step scaling func- 
tion to the continuum. 



Obviously there is very little cutoff" dependence, 
which is probably due to our use of a 1-loop 0(o) 
improved Wilson action. The extrapolation to 
the continuum, which on theoretical grounds is 
expected to be roughly linear in a/L for our back- 
ground field, presents no problem. Iterating the 
extrapolated values we gain information about 
a{s, 1.243) with s f« 2, 4, 8, 16, 24 [In the last step 
we used a factor 3/2 instead of 2, and the precise 
s values diff"er slightly from these integers as there 
are small mismatches in the iteration, which are 
easily taken into account. Moreover, all numbers 
have statistical errors. We omit these details in 
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this text, they can be found in [7] and are in- 
cluded in all figures]. The largest coupling oc- 
curring is g^{smaxL) = 3.48 with Smax ~ 24 and 
f{L) = 1.24. We expect 

Smax L — Lraax to be in 

the range of nonperturbative scales. Its precise 
relation to such a scale has to be determined in 
the last step. 

3. PHYSICAL SCALE AND COUPLING 

Rather than setting the scale through the string 
tension, we use the recently proposed [8] unit of 
length ro defined from the interquark force by 

rlF{ro) = 1.65. (2) 

If the force is identified with the one used in char- 
monium potential calculations, then ro ~ 0.5 fm 
is suggested. Advantages of ro over the string 
tension are discussed in [8]. On the basis of po- 
tential data from the literature, Lmaxlfo assumes 
the value 0.674(50) after extrapolating o ^ [7]. 
Our results are summarized in Fig. 2. The pertur- 
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Figure 2. Values of the running coupling together 
with perturbative evolution to 1-loop (dotted) 
and 2-loop (dashed), and a fit (solid). 

bative evolutions start from the smallest coupling 
and the fit is made with an effective 3-loop term. 
The error in the overall scale, which corresponds 
to a horizontal shift of all points and curves in the 
semi-logarithmic plot, is not included in Fig. 2. 



At this point we have learned that g'^{L = 

0. 027.3)ro) = 1.243. It remains to pass from our 
finite volume coupling to which by convention 
refers to the MS-scheme. At g = we should 
be well in the perturbative regime and can use 
the 1-loop formula [7] 

«Ms(9) = + 1-2556 a{qf + 0{a{qf). (3) 

The final result at the highest energy is now 

«Ms(9) = 0.1108(23)(10) at g = 37 (4) 

which corresponds to about 15 GeV. The second 
error is the estimated effect of missing higher or- 
ders in (3) and should be eventually eliminated 
by a 2-loop calculation, while the first one com- 
bines all other errors. As is obvious from Fig. 2, 
apart from computing (4) we have demonstrated 
that perturbative formulas should be trustworthy 
to evolve a^q) to yet higher energies. The desired 
link between low and high energy has thus been 
constructed. We plan to extend the method to 
dynamical fermions in the future. 
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